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MARTIAN ATMOSPHERIC LEE WAVES 

Joseph A. Pirraglia 
Laboratory for Planetary Atmospheres 
Goddard Space Flight Center 
Greenbelt, Maryland 20770 

ABSTRACT 

The Mariner 9 television pictures of Mars showed areas of extensive moun- 
tain lee wave phenomenon in the northern mid-latitudes during winter, hi most 
cases the characteristic wave length of the lee waves is readily observable and 
in a few cases the boundaries of the wave patterns, as well as the wave length, 
observed. The cloud patterns resulting from the waves generated by the flow 
across a mountain or crater are dependent upon the velocity profile of the air 
stream and the vertical stability of the atmosphere. Using the stability as inferred 
by the temperature structure obtained from the infrared spectrometer data, a two 
layer velocity model of the air stream is used in calculations based on the theory 
of mountain lee waves. The parameters that yield a pattern similar to that in a 
picture with a well defined wave configuration are a lower 3.1 km deep air stream 
of 40m/sec and an upper air stream of 85 m/sec. This calculation appears to be an 
upper limit of the wind speeds, with most of the pictures implying wind speeds of 
the lower layer to be less than 40 m/sec. These results yield magnitudes gener- 
ally in agreement with circulation models, in particular, the Leovy and Mintz 
(1969, J. Atmos. Sci. , .26, 1167) two layer numerical model. Under not too dif- 
ferent conditions they calculate winds of approximately 30 and 70 m/sec in the 
lower and upper layers respectively. 


MARTIAN ATMOSPHERIC LEE WAVES 


INTRODUCTION 

Atmospheric waves in the lee of mountains are an often observed and fairly 
well understood phenomena. The theory of mountain lee waves has been dis- 
cussed in a series of papers by Scorer (1949, 1953, 1954), among others, and 
in particular for an isolated peak by Scorer (1956) and Scorer and Wilkinson 
(1956). Observation of this familiar terrestrial phenomenon on another planet, 
as was the case during the Mariner 9 mission, greatly aids in the interpretation 
of the data. 

One of the objectives of the Mariner 9 mission to Mars was the investiga- 
tion of the atmospheric thermal and dynamical behavior. The thermal structure 
was obtained for extensive, though not complete, coverage in local time and 
latitude in the earlier part of the mission and for somewhat less extensive cov- 
erage in the latter part of the mission (Hanoi, et al. , 1972; Conrath, et al. , 
1973). The dynamical nature of the atmosphere can be investigated using tho 
thermal structure or can be inferred from apparent cloud motions (Briggs and 
Leovy, 1974). However, when the temperature data is supplemented by an 
additional observable, in particular, cloud patterns as in lee waves, winds can 
be inferred within the constraints imposed by the problem and serve as bounds 
on independently derived motion. This approach is discussed briefly by Briggs 
and Leovy (1974) and is the problem which will be addressed here. 
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During revolutions 179 to 220 of Mars the spacecraft was able to view the 
northern mid-latitudes during late winter. The pictures obtained during that 
period showed extensive waves in the clouds that comprise the edgo of the polar 
hood (see Briggs and Leovy for a collection of cloud pictures), hi many instances 
the cloud formations were obviously associated with atmospheric waves in the lee 
of mountains or craters. While these waves may be of no great consequence in 
the global circulation thoy do serve as a diagnostic means of determining wind 
velocities. We shall consider two examples of wave patterns and attempt to 
determine the air streams which, when used in calculations based on the theory 
of mountain lee waves with the vertical stability as inferred by the thermal 
structure, yield wave patterns similar to those observed. We shall sketch the 
development of Scorer, repeating parts with minor modifications for continuity. 

A simple extension of Scorer’s analysis will be used to introduce effects due to 
the horizontal dimensions of the crater generating the wave patterns. Rather 
than trying to reproduce die cloud patterns by calculating the vertical motions 
and condensation we shall calculate only the surfaces of constant phase. The 
wave patterns will be inferred by considering sections along the wave crest 
where the wave amplitude exceeds some reasonable but arbitrary prescribed 
value. 

hi most of the pictures the only easily observable parameter of the wave 
pattern is the wave length. While the whole pattern is needed to define the 
airstream profile the wave length alone does permit the establishment of 


limits on the alrstream near the surface. This will become apparent in the 
course of our exposition. 

It will be shown that there is agreement between the observations and 
theory and that in the latitudes, local time and season of the observations the 
theory does not suggest supersonic winds as might be inferred by the cloud 
patterns observed in some cases. 

FREE MODES OF STATIONARY WAVES 

The lee wave pattern is due to the discrete spectrum of waves generated 
by a crater ridge in an airstream. The discrete spectrum is the set of 
eigenmodes whose wavenumbers are the solutions to the zero frequency dis- 
persion equation which is derived from the homogeneous equation for the air- 
stream. 

The free modes arc infinite in lateral extent, beingthe type of wave generated 
by an infinitely long ridge not necessarily normal to the airstream. A lee wave 
pattern of a crater is obtained by summing over a parallel set of ridges oriented 
at all angles between ± n/Z from the normal to the airstream. Relative ampli- 
tudes of the modes are obtained by expanding a ridge in terms of the stationary 
waves with the ridge acting as a forcing function. The method is identical to 
a Green's function approach in which the ridge is used as the source term. 

The velocity, density, entropy and acoustic velocity of an undisturbed air- 
stream are designated by u 0 ,p 0 ,i? 0 and c Q respectively. Using the approach of 
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Eckart (19G0) the parameters associated with the kinematics and dynamics of 
the atmosphere are defined by 


(a, r .-L^ + J_ 

2 pQ C Q C 0 T 


<« N2 = t I ■* = 8 ii + r?) 3 


V 


( 1 ) 


„ r 
"P 


where the prime indicates differentiation with respect to z. y and R are the 
ratio of the specific heats and the gas constant, N is the Brunt- Vais ala fre- 
quency, r is a reciprocal scale height related to the adiabatic changes in the 
vertical direction and f3 is the static stability. 

With the stationary wave vector and horizontal space vector defined by k 
and r the vertical velocity, w, of a stationary wave is defined by 


w = W(p 0 c 0 )'^ exp 


lk ^rf 


z 

rdx 


( 2 ) 


and assuming the acoustic velocity to be greater than the stream velocity, the 
homogeneous equation for W is 


' N 2 | k |2 k. u„ k - ul n2\ 
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To obtain (3) it has also been assumed that the Rossby number is large and that 
the perturbation velocity of the stream is small compared to the stream velocity. 


One of the simplest models, and the one that we shall use, is a two layer 
model in which the velocity is uniform within a layer and in the same direction 


in both layers. The density decreases exponentially with altitude, the Brunt- 
Viiisala frequency is constant and both are continuous across the interface of 
the two layers. In addition, for the thermal structure of the atmosphorc in 
which we oi' e interested, the fifth and sixth terms aro much smaller than the 
first term in the brackots of (3) and may be neglected. Using the subscripts 
1 and 2 to identify quantities in the lower and upper layers respectively, taldng 
the interface at z = 0, the surface at z = -h and lotting 0 represent the angle 
between the wave vector and stream direction we have from (3) 

(a) W', 
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The interface conditions are 
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and at the surface the boundary condition is 

W, (-h) * 0. < 6 ) 

Solving the set (4), (5) and (0) under the condition that the wave energy 
decreases with height, that is the waves are trapped (Corby and Sawyer, 1958; 
Scorer, 1949), the eigenvalue equation for the wave numbers k is obtained, 
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The Scorer parameters lj and L are the limits of the possible values of k. 
Solutions of {7) give the valuos of the wavenumber k r of the free modes of the 
airstream. 

The roots k r define the modes of the stationary waves with wave vectors 
at an angle <p to the direction of tho airstream. The modes have the horizontal 
behavior, 


w ~ e ik f (xcos0 + y sin <p) (9) 

where x is in the direction of the airstream and y is normal to the airstream. 

If the waves are assumed to be generated by a ridge at an angle <p to the stream 
and at a distance d from the origin of the x-y plane, as shown to Fig. 1, we 
have, 

w ~ c ik r licos(0-i£) + d] =e*M (10) 


with r and P as shown to Fig. 1. 
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Solutions of the form (10) represent the free modes but in order to deter- 
mine the amplitudes the source of the wtives must bo considered. In addition, 
wo are interested in the vertical excursion of an atmospheric parcel not its 
vertical velocity since the clouds will presumably form or bo moro dense 
along the crests of the waves in the streamlines. The vortical excursion, 
of a streamline for the wave of wavenumber K is related to the vertical 
velocity by 


(z ' {) ■ #{ < u > 

to within a constant phase parameter. The ridge is described in terms of the 
stationary waves by 

£ ($) = a (l + — ^ = — Real pari f c* kb + ik * dk. ■ 

* \ b 2 / n J o 


Normalizing to the amplitude of the displacement at the surface and using 
tliis with (112) which represents the boundary condition at z = -li when a ridge 
is in the ailrstream we have, with He designating the real part, 
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represents the vertical displacement of a streamlJr-j duo to the stationary 
waves with wave vectors porpendicular to the riuge and at an angle <p to tho 
airstroam as shown in Figure 1. 

If wo now toko 


(a) 

(b) 


£ w = rcos(0-0) + d 
£ L “ r cos (0 - 0) - d 

(14) 

, z 

A - ab lp 0 (z)/p 0 (-li)J exp (— f rdzj , 

(15) 


consider a set of parallel ridges separated by 2d, and integrate (13) over 0 

, 7T 7 r 

from - ’T" to ~ we have 

a 4 


ffer , 0 ) = - 
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/l0 /dk~|^ o'^Ce^w+e^L) (16) 


where f (z, r, 0 ) represents the complete solution for the displacement of the 
streamlines of flow across a crater. The discrete free modes represented by 
the solutions to (7) are at the poles of the integrand. Integration in the complex 
k-plane is chosen so that the waves are downstream from the crater, the inte- 
gration paths being different depending upon whether £ w or £ L are less than or 
greater than zero, (See Scorer, 1949, for a mo.:e detailed discussion of the 
k-plane integration.) Evaluation of the lee waves is straight-forward and is 
simply the residue of the poles; 1 
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whore j * 1 when -h < z < 0 and j = 2 when z > 0. 


(17) 


Because of the e' k r b term tlie mode with the smallest root will bo dominant and 


we shall consider it only and drop the summation and the r-subsoript. 

Separating the intogral into parts containing £ w and £ L the integral is eval- 
uated by the method of stationary phase where the stationary points <j>„ and 4> L 
are defined by 

(a) j (kt.) - 0 

(18) 

(b) ± <ky ■ 0 


On constant phase surfaces defined by the constant M, 

(a) kfc w = M (19) 

(b) k£ L = M . 

Using (18) and (19) equations are obtained for 0 and r in terms of the stationary 
points <f> w and 0 L ; 
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where (1/k) 9k/90 is obtained from the dispersion equation (7) and the upper 
and lower signs are associated with the W and L subscripts respectively. 
Equations (20) and (21) yield r and along constant phase fonts where the 
phrase is determined by the constant M. 

The stationary phase evaluation of (17) gives the result 
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where all the terms are evaluated at the stationary point 0 s , which equals 0 or 

depending upon whether the wave is due to the windward or leeward rim of 

3 2 

the crater and the sign of tt/ 4 is + or - depending upon whether (k£ , ) 

30 2 W ’ 

is greater than or less than zero. When the second derivative is equal to zero 
(22) is singular and the solution is similar to (22) but with the curly bracketed 
term replaced by 


np (±\ 
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with T (4/3) being the gamma function. 
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If for 3 2 (Ic£ w l )/3i p 2 > 0 wc take M= (G + 8n -1) and for a 2 k k^ w>L )/30 2 < 0 

take M «= {6 + 8n + 1) ~ with n - 0, 1, 2, , the constant phase surfaces are 

alternately the wave crests duo to the windward and leeward rims of the crater. 
Using (20), (21) and (22) the amplitude and coordinates of the wave crests can be 
calculated as a function of <f> s . 

DISCUSSION OF OBSERVATIONS AND CALCULATIONS 

The vertical displacement of the stream lines in addition to being dependent 
upon the stream parameters also depend upon the parameters a and b which 
describe the crater ridge, as is evident in (15) and (22). As pointed out pre- 
viously, to determine the wave pattern it is not necessary to calculate the wave 
amplitudes over the total wave field but all that is needed is the contours of the 
wavo crests since the clouds will form initially and be most dense along a crest. 
If the height and slope of the ridge is assumed to be small enough to satisfy the 
condition that the perturbation velocity is small as compared to the stream 
velocity then the horizontal contours of the wave crests are independent of the 
vertical scales of the crater ridge leaving only the easily observable crater 
diameter, in addition to the airstream parameters, needed to determine a wave 
crest pattern. However, the width of the ridge as compared to the wave number 
does have an influence on the observed details of the cloud pattern. The factor 
e' kb in (22) depends upon the angle between the airstream and the normal to the 
crater at its point of maximum contribution to the wave. This implies that the 


l'idgo width determines, to some extent, the sections along the crest where the 
amplitude exceeds that required to cause condensation of the clouds. Thus, 
the location of the clouds along a wave crest is somewhat dependent upon the 
ridge width. 

Figure 2 shows a Mariner 9 television picture taken during the 185th revolu 
tion of the planet by the spacecraft. The 100 lan dia. frost coated crater in the 
lower left comer is at 55N latitude and 145W longitude, and the evening termi- 
nator runs across the upper left comer of the picture. At the time of the pic- 
ture the sub solar latitude was 7°S, making it late winter in the northern hemi- 
sphere; The clouds in the picture are composed of water ice and are along the 
edge of the north polar hood which was receding at that time. (The discovery 
of H 2 0 ice clouds in the Martian atmosphere is discussed in Curran, et al. , 
1973. ) The apparent wind direction is along the dark streak extending east of 
the crater. Two sets of overlapping waves can be seen in the southern part of 
the wave system generated by the crater whom it is relatively free of inter- 
ference from other waves. There is a transverse set normal to the wind 
direction directly in the lee of the crater and a div erging set with crests that 
intersect directly behind the crater. The waves fill a wedge shaped region 
whose edge makes an angle of approximately 18° with the apparent wind direc- 
tion and the wave length of the transverse waves along the centerline of the 
pattern is approximately 60 km. 


Using the model described in the previous section by assigning values to 
N \ which depends upon the temperature profile, to the airstream speeds u^ 
and u 02 and to the depth h of the lower airstream we attempt to construct the 
wave pattern which could generate the cloud configuration of Figuro 2. Unfor- 
tunately, there is no temperature data available for the time of the picture but 
other data for the same latitude and local time Indicates that the atmosphere 
was substantially subadiabatic and in fact may have had temperature inversions. 
As an approximate mean, in our models we shall use an isothermal atmosphere 
where from (lb), (3 - j3 iso ^ g/(c p T). Any stream velocities can then be considered 
to be normalised to the isothermal case and the velocity associated with any 
other value of (3 is obtained by multiplying the isothermal case by (|3/j3 iso ) Vl . 

Figure 3 shows families of airstream models as a function of wedge-angle 
and transverse wavelength. Part (a) of the figure is for lower layer velocities 
of 20 and 40 m/s and part (b) is for 30 and 50 m/s. In each part of the figure 
lines of constant lower layer depth and upper layer velocity are plotted. The 
point associated with the pattern of Figure 2 is indicated in both parts. The 
dotted portion of each section is the limit beyond which there is no transverse 
wave across the center line of the pattern. Evidently, any one of an infinite but 
bounded set of airstreams produce the correct angle and wavelength. An obvious 
one from the figure is the case with u 01 - 40 m/s, u 02 = 85 m/s and h = 11 Ion. 
Any model with u 01 between approximately 30 and 45 m/s and with the appropri- 
ate values of u 02 and h could satisfy the two constraints. This does not imply 
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a laok of uniqueness but rather that the whole contour of the wave must be 
matched, not only the two pieces of information contained in the wedge angle 
and wav ength along the centerline of the pattern. With our three parameter 
model we cannot expect to exactly duplicate the wave crest contour but simply 
choose that one wliich best fits the observation. 

Figure 4 shows three wave patterns calculated using the method outlined 
in the previous section. All have identical wedge angles and wavelengths along 
the centerline of the pattern. The cases five, with [u 01 , u 02 , h] representing 
the three parameters, [35, 83, 8 J , [ 40, 85, 11] and [45, 92, 15] as indicated. 
These were chosen by interpolation in Figure 3. Model [35, 83, 8] has wave 
crests with more curvature than those in Figure 2 and model [45, 92, 15 ] are 
too flat while [ 40, 85, 11] is a reasonably good reproduction. The observed 
patterns of Figure 2 lie between the two extreme models of Figure 4 which 
yield velocities without very large differences so that the gross features of the 
airstream would be obtained in any case. 

The clouds would form along the contours of Figure 4 depending upon the 
shape of the ridge. Those in Figure 2 appear to be along the transverse waves 
of Figure 4 associated with the outer wedge, which is due to the windward ridge, 
and along the diverging waves of the inner wedge, which is due to the lee ridge. 
This suggests that the windward ridge of the crater is not as steep as the lee- 
ward ridge. 
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Figure 5 is an example of a less common wave pattern. Tills picture covers 
approximately the same area as Figure 2 and the cloud pattern is due to the same 
crater but proceeds that of Figure 2 by one Martian day. It is apparent that 
there is no transverse wave system. This is because waves cannot form in that 
direction for the conditions of the airstrenm at that time since, as predicted by 
(7), there is no solution for lc when absolute values of <p are less than some non- 
zero magnitude. As before there is more than one model that will give the 
proper wedge angle and wavelength along some prescribed direction. In tills 
c ase the angle is readily apparent but the wavelength is not and the best we can 
do is to match the angle and roughly approximate the wavelength along the 
streaks. A candidate airstreara is [20, 85, 4 1 whose lee wave pattern is shown 
in Figure 6, There is a good similarity between Figures 5 and 6 especially in 
the first three waves in the lower part of the wake and in the characteristic 
streaks flaring out from the crater. The air stream of Figure G followed by that 
of Figure 4 by one day implies an upper layer of 85 m/s, uniform in time, over 
a lower layer of 20 m/s velocity increasing to 40 m/s, and of 4 km depth increas- 
ing to 11 lan. This behavior is plausible being suggestive of lower layer growth. 

The theory is essentially good only for far fields, i.e. , icr large, but can 
in fact be used to give near field results for the free modes while ignoring the 
continuous part of the wave spectrum. An example of the agreement between 
theory and observation is indicated by Figure 7. This picture shows the first 
and second wave crests of the lee waves of a 45 km crater at 63N, 347W (this 
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picture is 200° in longitude from, and 3-^ daya prior to Figure 2). The theory 
predicts that the first crest is 5/8ths of a wavelength from the windward ridge 
of the crater. The wavelength between the two crests is 45 1cm and the distance 
from the windward rim to the first crest is 28 Ion which is in very good agree- 
ment with the theory. 

Returning to Figure 3 for a moment and confining the discussion to patterns 
with complete transverse waves, the changing character of the wave pattern 
can be observed as a function of the airstream parameters. As the Imvnr 
stream velocity increases the minimum depth of the layor must increase to 
produce patterns with complete transverse waves. For increasing lower stream 
depth the pattern becomes less dependent upon the upper stream velocity as, for 
example, alongtheh = 10 line for u Q1 = 20. As the difference between the upper 
and lower stream velocities decreases the pattern becomes less dependent upon 
the depth as is seen along the u 02 = 30 line for u Q1 = 20, Both cases are intuitively 
evident. Of particular interest to us is the observation that for a given lower 
layer velocity a lower limit can be placed on the wavelength by the constraint 
that k cannot exceed 1 lf as seen from (8). k equals lj where the contours of 
Figure 3 touch the zero wedge angle line. The limit on k is established without 
the aid of Figure 3 but the figure shows that for a given wavelength wide wedge 
angles infer small lower stream velocities. The implications of which in con- 
junction with the observations will be discussed subsequently. 
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NONUNIFORM STREAM DIRECTIONS 

If the two streams of the model have different directions we have another 
parameter of the airstream that must be determined from the observations. 
Because, for the most part, wo cannot even determine three parameters no 
attempt will be made to determine a fourth as expressed by the angle between 
the upper and lower streams. However, calculations shoving the types of 
patterns that are the result of nonuniform stream directions will be presented. 
If the angle between the two streams is p tie eigenvalue equation for k is, 


v cot Oh)-i 
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u 02 

with (23) and (24) used in the place of (7) and (8) the analysis proceeds as befox'e. 
Using a model with a lower stream velocity of 20 m/s and depth of 5 lan and 
an upper stream velocity of 80 m/s the patterns were calculated for cases in 
which the angle between the streams was 45° and 90°. Similar calculations 
were made with a 10 Ian deep lower stream of 40 m/s. The results are shown 
in Figure 8. The calculations indicate that a pair of streams with similar depth 
and wind speeds will yield shorter wavelengths for larger angles between the 
streams. However, the observed wave patterns generally are not suggestive of 


/s \\ 
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the type of patterns shown In Figure 8 and the conclusions drawn from the parallel 
stream case of the previous section should be valid. The streams may bo off 
parallel but not by the amounts used for the calculations of Figure 8. If there 
were indications of greatly different wind directions, without observations of 
isolated lee wave patterns it would be exceedingly difficult to determine four 
parameters of the airstream. 

CONCLUSIONS 

Regardless of the simple model used for the airstream the calculations 
appear to be capable of reproducing the observed wave patterns. Insofar as 
the model is a reasonable approximation to the actual airstream, were it not 
for the paucity of temperature data at the latitude and season of the observations, 
accurate assessments of the wind velocities could be made. Since all our cal- 
culations were for a 2C0K isothermal atmosphere the velocities would have to be 
adjusted by the square root of the ratio of the stabilities, as was pointed out 
previously. If the temperature profile approached an adiabatic state the velocities 
as calculated would decrease toward zero. If there were a temperature inver- 
sion, the velocities would increase, but to increase the velocities by 25% over 
those of the isothermal case would require an inversion of more than GI</Km. 

Such an inversion is large and not very likely, making the possibility of velocities 
greater by 25% improbable. 

The first wave in the clouds shown in Figure 2 is not normal to the apparent 
stream direction and may be due to the nonuniformity of the wind direction with 


altitude. As illustrated previously, the wave patterns are not symmetric about 
a horizontal lino when tho windstreams have different directions but the wavo 
patterns do not suggest a large difference in the wind directions. With small 
differences and without isolated patterns it is not possible to determine the dif- 
ference in the stream direction but^noverthclcss, one may use the wavelengths 
to give an estimate of the winds in the lower layer. 

In most of the televh ion pictures the wavelengths of the apparent Ice waves 
are approximately 30 km and of the type shown in Figures 2 and 4 rather than 
those of Figures 5 and 6. Although the lack of many clear cases of wave patterns 
showing the wedge angle prevents a determination of the upper airstream velocity, 
an upper limit can be put on the velocity of the lower stream. From Figure 3 
it is apparent that wavelengths of 30 km are associated only with airstreams with 
lower stream velocities loss than 40 m/s since the observations suggest wide 
wedge angles and from the discussion in the previous section they are associated 
with the smaller lower stream velocities for a given wavelength. 

Pictures like Figure 5 give the impression of very high velocities but the 
results of the calculations shown in Figure 6 and compared with Figure 5 imply 
that there is no need to invoke extremely high velocities to explain the particu- 
lar wake pattern and, in fact, the surface velocity is less than the more common 
cases with transverse waves. 

The conclusion may be made that the many observations of lee waves sug- 
gest that the near surface wind speeds, except possibly on a very localized 
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scale, do not greatly exceed, and most of the time are loss than, 40 m/scc in 
the northern midlatitudos during the late winter. This result is in - 'roement 
with the wind magnitudes obtained independently, in particular with the two 
level numerical model of Leovy and Mintz (19G9) where they calculated the 
mean zonal winds at SON and at winter solstice to bo 30 m/s at 3Km and 70 m/s 
at 13. 5 Km. Winds of this magnitude are in substantial agreement with the 
winds implied by the lee waves. 

One may question the validity of the two layer model used to describe the 
wind profile. The two layer model with a discontinuous velocity profile is sub- 
ject to Kelvin-Helmholtz instabilities at all wavelengths. However, a three 
layer model with uniform upper and lower velocities connected continuously 
through the varying middle layer is unstable to only a small band of wave- 
numbers which are dependent upon the Richardson number and any Kelvin- 
Helmholtz instability wavelength is likely to be far removed from the lee wave 
lengths. If the product of the transition layer thickness and the wavenumber is 
much less than unity the two layer approximation in which the transition layer 
is absent is a good approximation for the calculation of lee waves in a contin- 
uously varying velocity profile airstream. The only question that remains is 
whether or not the three layer model is an adequate representation of the air- 
stream. The extensive periodic train of lee waves observed is the result of 
the waves being trapped in a waveguide formed by layers in which the Scorer 
parameter, 1 = N sec0/u Q , allows propagation in a lower layer and reflection 


in an upper layer. The throe layer model is goneral enough to permit trapping 
of the wavos and continuity of the velocity profile and should tlion bo at least a 
rough approximation to the actual airstream and its two layer approximation 
is adequate for the calculation of lee waves. 
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FIGURE CAPTIONS 

Figure 1. Horizontal geometry of airstream, ridge and wave directions. The 
ridge is at a distance d from the x-y and r -0 coordinate origin and 
its normal is at an angle 0 to the airstream u 0 . The wave vector k 
is normal to the ridge. 

Figure 2. Mariner 9 television photograph taken during the 185th revolution 
of the planet by the spacecraft showing an extensive system of lee 
waves. The white circular region in the lower left comer is a 
frost coated 100 Ion diameter crater located at 55N latitude and 
145W longitude. The evening terminator runs diagonally up from 
the lower right comer. 

Figure 3. Plots of wedge angle vs. wavelength of transverse wave for various 
upper and lower airstream speeds and lower layer depths. Part 
(a) shows the two surfaces with lines of constant depth and upper 
stream speeds for lower stream speeds of 20 and 40 meters per 
sec. Part (b) shows similar surfaces for lower stream speeds of 
30 and 50 meters per sec. The wavelength and wedge angle of the 
wave system of Figure 2 is indicated by the black circle in both 


parts of the figure 
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Figure 4. Calculated wave patterns. Part (a) is the wave system generated by 
a 100 km dia, crater in an airstream with a 45 m/s, 15 km deep 
lower stream and a 02 m/s upper stream. Part (b) is with a 40 m/s, 
11 1cm deep lower stream and 85 m/s upper stream. Part (c) is 
with a 35 m/s, 7. 8 km deep lower stream and 82. 5 m/s upper 
stream. 


Figure 5. Mariner 9 television photograph taken during revolution 18,3, The 
crater on the loft border is the same one shown in Figure 2. The 
time of this picture precedes that of Figure 2 by two Martian days 
and there are no transverse waves here as are evident in Figure 2. 


Figure G. The calculated wave system generated by a lower stream of 20 m/s 
and 3. 5 km depth and an upper stream of 85 m/s blowing over a 
100 1cm dia. crater. 

Figure 7. Mariner 9 picture of two lee wave crests due the 45 km dia. crater 
at 63N latitude and 347W longitude. 

Figure 8. Calculated wave patterns for two layer air streams in which the 
upper and lower velocities have different directions. The speeds, 
lower layer depth and angle between the velocities of the upper and 
lower streams are indicated in each section of the figure. 
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Figure 2. Mariner 9 Television Photograph Taken During the 185th Revolution 
of the Planet by the Spacecraft Showing ;ui Extensive System of Ia*o 
W aves. The White Circular Region in the Lower Loft Corner is a 
Frost Coated lot) km Diameter Crater Located at 55N Latitude and 
145W Longitude. The Evening Terminator Runs Diagonally lip 
From the Lower Right Corner. 
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Figure 3. Plots of wedge angle vs. wavelength of transverse wave for 
various upper and lower airstream speeds and lower layer 
depths. Part (a) shows the two surfaces with lines of constant 
depth and upper stream speeds for lower stream speeds of 20 and 
40 meters per sec. Part (b) shows similar surfaces for lower 
stream speeds of 30 and 50 meters per sec. The wavelength and 
wedge angle of the wave system of Figure 2 is Indicated by the 
black circle in both parts of the figure. 
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Figure 4. Calculated wave patterns. Part (a) is the wave system generated 
by a 100 km dia. crater in an airstroam with a 45 m/s, 15 km 
deep lower stream and a 92 m/s upper stream. I*art (b) is with 
a 40 m/s, 11 km deep lower stream and 85 m/s upper stream. 
Part (c) is with a 35 m/s, 7. 8 km deep lower stream and 82. 5 
m/s upper stream. 



Mariner I) Television Photograph Taken During Revolution lHM. 
The Crater on the Left Border is the Same Cine Shown in Figure 2 
The Time of this Pieture Precedes that of Figure 2 by Two 
Martian Days and There are No Transverse Waves Here as are 
Kvident in figure 2. 
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Figure 8. Calculated wave patterns for two layer air s .reams in which the 
upper and lower velocities have different df lections. The speeds 
lower layer depth and angle between the .'flocities oi the upper 
and lower streams ~ ,T >dlcated In each secHou of the figure. 



